International Journal of Modern Physics C, Vol. 0, No. (2000) 000-000 
© World Scientific Publishing Company 



Dynamical aspects of equilibrium and off-equilibrium simulations 
of a three-dimensional Heisenberg spin glass 



H. Kawamura 

Kyoto Institute of Technology, Sakyo-ku, Kyoto 606-8585, Japan 



and 



K. Hukushima 

ISSP, University of Tokyo, Minato-ku, Tokyo 106-8666, Japan 



Spin-glass and chiral-glass orderings of a three-dimensional isotropic Heisenberg spin 
glass are studied both by equilibrium and off-equilibrium Monte Carlo simulations with 
emphasis on their dynamical aspects. The model is found to exhibit a finite-temperature 
chiral-glass transition without the conventional spin-glass order. Although chirality is 
an Ising-likc quantity from symmetry, universality class of the chiral-glass transition ap- 
pears to be different from that of the standard Ising spin glass. In the off-equilibrium 
simulation, while the spin autocorrelation exhibits only an interrupted aging, the chiral- 
ity autocorrelation persists to exhibit a pronounced aging effect reminiscent of the one 
observered in the mean-field model. 

1. Introduction 

Ordering of complex systems has attracted interest of researchers working in the field 
of numerical simulations. Well-known examples of such complex systems may be a 
variety of glassy systems including window glasses, orientational glasses of molec- 
ular crystals, vortex glasses in superconductors and spin-glass magnets. Often, in 
the dynamics of such complex systems, characteristic slow relaxation is known to 
occur. It has been a great challenge for researchers in the field to clarify the nature 
and the origin of these slow dynamics, as well as to get fully equilibrium properties 
by overcoming the slow relaxation. In particular, spin glasses are the most exten- 
sively studied typical model system, ipr which numerous analytical, numerical and 
experimental works have been madeBl 

Studies on spin glasses now have more than twenty years of history. Main focus 
of earlier studies was put on obtaining the equilibrium properties of spin glasses. 
While extensive experimental studies have clarified that the spin-glass magnets ex- 
hibit an equilibrium phase transition at a finite temperaturejd the true nature of 
the experimentally observed spin-glass transition and that of the low-temperature 
spin-glass phase still remain open problems. It has been known that the magnetic 
interactions in many of real spin-glass materials are Heisenberg-like, in the sense 
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that the magnetic anisotropy is much weaker than the exchange energy. In apparent 

opic Heisen- 
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Apparently, 



contrast to experiments, numerical simulations have indicated 
berg spin glass exhibits only a zero-temperature transitions 
there is a puzzle here: How can one reconcile the absence of the spin-glass order in 
an isotropic Heisenberg spin glass with the experimental observation? 

In order to solve this apparent puzzle, a chirality mechanism of experimentally 
observed spin-glass transitions was recently proposed by one of the authorspu on the 
assumption that an isotropic 3D Heisenberg spin glass exhibited a finite-temperature 
chiral-glass transition without the conventional spin-glass order, in which only spin- 
reflection symmetry was broken with preserving spin-rotation symmetry. "Chiral- 
ity" is an Ising-like multispin variable representing the sense or the handedness of 
the noncoplanar spin structures. It was argued that, in real spin-glass magnets, the 
spin and the chirality were "mixed" due to the weak magnetic anisotropy and the 
chiral-glass transition was then "revealed" via anomaly in experimentally accessi- 
ble quantities. Meanwhile, theoretical question whether there really occurs such 
finite-temperature chiral-glass transition in an isotropic 3D Heisenberg spin glassy 
crucial assumption of the chirality mechanism, remains somewhat inconclusive .Effi 

More recently, there arose a growing interest both theoretically and experi- 
mentally in the off- equilibrium dynamical properties_of spin glasses. In particu- 
lar, aging nh en o m en a. observed in many spin glassest^l have attracted attention of 
researchers JlJil3 Unlike systems in thermal equilibrium, relaxation of physical quan- 
tities depends not only on the observation time t but also on the waiting time t Wl 
i.e., how long one waits at a given state before the measurements. Recent studies 
have revealed that the off-equilibriu m dy namics in the spin-glass state generally 
has two characteristic time regimes. One is a short-time regime, to -C t <C t w 
(to is a microscopic time scale) , called "quasi-equilibrium regime" , and the other is 
a long-time regime, t 3> t w , called "aging regime" or "out-of-equilibrium regime". 
In the quasi-equilibrium regime, the relaxation is stationary and the fluctuation- 
dissipation theorem (FDT) holds. On the other hand, in the aging regime, the 
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I succeeded in reproducing 



relaxation becomes non-stationary, FDT broke 
off-equilibrium dynamics on Ising-like modelsc 
some of the features of experimental results, it is clearly desirable to study the dy- 
namical properties of a Heisenberg spin-glass model to make a direct link between 
theory and experiments. 

In the present article, we report on our recent results of equilibrium as well as 
off-equilibrium Monte Carlo (MC) simulations on an isotropic 3D Heisenberg spin 
glass, with emphasis on their dynamical aspects. Ordering properties of both the 
spin and the chirality are studied, aimed at testing the validity of the proposed chi- 
rality scenario of spin-glass transitions. We note that MC simulation is particularly 
suited to this purpose, since, at the moment, chirality itself is not directly measur- 
able experimentally. By contrast, it is straightforward to measure the chirality by 
numerical simulations. 



Dynamical aspects of a three-dimensional Heisenberg spin glass 3 



2. Model 

The model we simulate is the classical Heisenberg model on a simple cubic lattice 
with the nearest-neighbor random Gaussian couplings, Jy, defined by the Hamil- 
tonian 

n Y^.i.^.-s,. (2.i) 

<ij> 

where = (Sf , Sf, Sf) is a three-component unit vector, and the sum runs over 
all nearest-neighbor pairs with N = L x L x L spins. Jy is the isotropic exchange 
coupling with zero mean and variance J. 

Frustration in vector spin systems often causes noncollinear or noncopianar 
spin structures. Such noncollinear or noncopianar orderings give rise to_a nontrivial 
chirality according as the spin structure is either right- or left-handed. o We define 
the local chirality at the i-th site and in the /^-th direction, i f° r three Heisenberg 
spins by, 

Xin = Si+g„ • (Si x Sj_e,J, (2.2) 

where e M (fi = x, y, z) denotes a unit lattice vector along the //-axis. The chirality is 
a pseudoscalar in the sense that it is invariant under global spin rotation but changes 
sign under global spin reflection. Chiral order is related to a possible breaking of 
the spin-reflection symmetry with preserving the spin-rotation symmetry. 



3. Equilibrium simulations 

First, we report on our equilibrium MC simulations of a fully isotropic 3D Heisen- 
berg spin glass defined by Eq. (2.1)Ja Monte Carlo simulationJs performed based 
on an exchange method developed by Hukushima and NemotoJi3 where whole con- 
figurations at two neighboring temperatures of the same sample are occasionally 
exchanged keeping the system at equilibrium. The method has turned out to be efH= 
cient in thermalizing the system showing slow relaxation over free energy barriers .Ej 
By this method, we succeed in equilibrating the system down to the temperature 
considerably lower than those attained in the previous simulations. We run in paral- 
lel two independent replicas with the same bond realization and compute an overlap 
between the chiral variables in the two replicas, 



5x = ^yX^XiJxi?- (3-1) 



The lattice sizes studied are L = 6,8,10,12,16 with periodic boundary con- 
ditions. In the case of L = 16, for example, we prepare 60 temperature points 
distributed in the range [0.10J, 0.25 J] for a given sample, and perform 1.2 x 10 6 
exchanges per temperature of the whole lattices combined with the same number of 
standard single-spin- flip heat-bath sweeps. For L = 16, we equilibrate the system 
down to the temperature T/J = 0.10, which is lower than the minimum tempera- 
ture attained previously. Sample average is taken over 1500 (L = 6), 1200 (L = 8), 
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Fig. 1. Log-log plots of the time dependence of the equilibrium spin (a) and chirality (b) 
autocorrelation functions at several temperatures. The lattice size is L = 16 averaged over 64 
samples. 



640 (L = 10), 296 (L = 12) and 136 (L = 16) independent bond realizations. 
Equilibration is checked by monitoring the stability of the results against at least 
three-times longer runs for a subset of samples. 

Information about the equilibrium dynamics can be obtained via the spin and 
chirality autocorrelation functions defined by 

= ^EKs^M-s^ + M}], (3.2) 

i 

= ^£[(x*(*o)Xi/*(* + to))], (3.3) 

where (■ • •) and [■ • ■] represent the thermal average and the average over bond disor- 
der, respectively, while t$ denotes some initial time at which the system is already 
fully equilibrated. Although we have used the exchange MC method for thermal- 
izing the system, the time evolution during the measurements of autocorrelations 
has been made via the standard single-spin-flip heat-bath updating. MC time de- 
pendence of the calculated C s (t) and C x {t) are shown in Figs, [l] on log- log plots 
for several temperatures. As can be seen from Fig. 0(a), C s (t) shows a downward 
curvature at all temperature studied, suggesting an exponential-like decay char- 
acteristic of the disordered phase, consistently with the absence of the standard 
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spin-glass order. In sharp contrast to this, C x {t) shows either a downward curva- 
ture characteristic of the disordered phase, or an upward curvature characteristic of 
the long-range ordered phase, depending on whether the temperature is higher or 
lower than T/J ~ 0.16, while just at T/J ~ 0.16 the linear behavior corresponding 
to the power-law decay is observed: See Fig. |l|(b). Hence, our dynamical data indi- 
cates that the chiral-glass order without the standard spin-glass order takes place 
at Tqq/J = 0.160 ± 0.005, below which a finite chiral Edwards- Anderson (EA) 
order parameter q^, > develops. From the slope of the data at T = Tqg, the 
exponent A characterizing the power-law decay of C x (t) t~ x is estimated to be 
A = 0.193 ± 0.005. We note that the occurrence of a chiral-glass transition was also 
supported by the behaviors of the static Binder ratios (not shown here) .! 
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Fig. 2. Chiral-overlap distribution function below Tqq. The temperature is T/J = 0.1. 




Fig. 3. MC-time evolution of the chiral-overlap distribution function at T/J = 0.1 with L = 16 
(left). MC-time dependence of the height of P(q x ) at g x = and at q x ~ qea (right). 



Establishing the existence of a finite-temperature chiral-glass transition, we pro- 
ceed to the study of the properties of the chiral-glass ordered state itself. In Fig. [|, 
we display the distribution function of the chiral-overlap defined by 



P(q , x ) = [(S(q x -q' x ))}, 



(3.4) 



6 H. Kawamura & K. Hukushima 



calculated by the exchange MC method at a temperature T/J = 0.1, well below 
the chiral-glass transition temperature. As is evident from Fig. the shape of 
the calculated P(q x ) is somewhat different from the one observed in the standard 
Ising-like models such as the 3D EA model or the mean-field SK model. As usual, 
P{q x ) has standard "side-peaks" corresponding to the EA order parameter ±<7cG' 
which grow and sharpen with increasing L. In addition to the side peaks, an 
unexpected "central peak" at q x = shows up for larger L, which also grows and 
sharpens with increasing L. This latter aspect, i.e., the existence of a central peak 
growing and sharpening with the system size, is a peculiar feature of the chiral- 
glass ordered phase never observed in the EA or the SK models. Since we do not 
find any sign of a first-order transition such as a discontinuity in the energy, the 
specific-heat nor the order parameter q^Q, this feature is likely to be related to a 
nontrivial structure in the phase space associated with the chirality. We note that 
this peculiar feature is reminiscent of the behavior characteristic of some^mean-ficld 
models showing the so-called one-step replica- symmetry breaking (RSB).tJ In Figs. ||, 
we show the (exchange) MC time dependence of P(q x ), together with those of the 
height of the central peak P(0) and of the side peak P(q , ea)- These quantities 
are found to reach stationary values exponentially fast with the correlation time of 
order of 10 5 MC steps, while actual spin configurations continue to change via the 
temperature-exchange process and the relatively rapid relaxation realized at higher 
temperatures. 

4. Off-equilibrium simulations 

In this section, we report on our results of off- equilibrium MC simulations .El Unlike 
the case of equilibrium simulations, the system here is never in full thermal equi- 
librium. Recent studies have revealed that one can still get useful information from 
such off-equilibrium simulations, even including certain equilibrium properties. The 
quantities we are mainly interested here are the off-equilibrium spin and chirality 
autocorrelation functions defined by 



MC simulation is performed based on the standard single spin-flip heat-bath 
method here. Starting from completely random initial configurations, the system 
is quenched to a working temperature. Total of about 3 x 10 5 MC steps per spin 
are generated in each run. Sample average is taken over 30-120 independent bond 
realizations, four independent runs being made using different spin initial conditions 
and different sequences of random numbers for each sample. The lattice size mainly 
studied is L = 16 with periodic boundary conditions, while in some cases lattices 
with L = 12 and 24 are also studied. 




(4.1) 




(4.2) 
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The spin and chirality autocorrelation functions at a low temperature T/J = 
0.05 are shown in Fig. |as a function of t. For larger t w , the curves of the spin 
autocorrelation function C s come on top of each other in the long-time regime, 
indicating that the stationary relaxation is recovered and aging is interrupted. This 
behavior has bee^fjxnected because the 3D Heisenberg spin glass has no standard 
spin-glass order .□'tnj'B'BLlla By contrast, the chiral autocorrelation function C x shows 
an entirely different behavior: Following the initial decay, it exhibits a clear plateau 
around t ~ t w and then drops sharply for t > t w . It also shows an eminent aging 
effect, namely, as one waits longer, the relaxation becomes slower and the plateau- 
like behavior at t ~ t w becomes more pronounced. It should be noticed that the 
plateau-like behavior observed here has been hardly noticeable in simulations of the 
3D EA model. 
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Fig. 4. Off-equilibrium spin (a) and chirality (b) autocorrelation functions at a temperature 
T I J = 0.05 plotted versus log 10 t for various waiting times t w . The lattice size is L = 16 averaged 
over 66 samples. 



While the plateau-like behavior observed in C x is already suggestive of a nonzero 
chiral Edwards-Anderson order parameter, q^,Q > 0, more quantitative analysis 
similar to the one recently done by Parisi et al for the 4D Ising spin glassEi is 
performed to extract q^Q from the data of C x in the quasi-equilibrium regime. 
Finiteness of is also visible in a log-log plot of C x versus t as shown in the inset 
of Fig. U, where the data show a clear upward curvature. We extract <7qq by fitting 
the data of C x for t w — 3 x 10 5 to the power-law form, C(t w , t + t w ) = q EA + in 
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Fig. 5. Temperature dependence of the Edwards- Anderson order parameter of the chirality. The 
data are averaged over 30-120 samples. Inset exhibits the log-log plot of the t-dependence of the 
chiral autocorrelation function in the quasi-equilibrium regime for L = 16 and l„ = 3x 10 5 . 

the time range 40 < t < 3,000 satisfying t/t w < 0.01. The obtained g^Gi plotted 
as a function of temperature in Fig. |5| clearly indicates the occurrence of a finite- 
temperature chiral-glass transition at Tcg/ J = 0.157 ± 0.01 with the associated 
order-parameter exponent flcG — 1.1 ± 0.1. The size dependence turns out to be 
rather small, although the mean values of q^Q tend to slightly increase around 
Tcg with increasing L. Since both finite-size effect and finite-tu, effect tend to 
underestimate Qqq, one may regard the present result as a rather strong evidence 
of the occurrence of a finite-temperature chiral-glass transition. 

The present estimate of the transition temperature Tcg/ J = 0.157 ± 0.01 is 
in very good agreement with the equilibrium estimate in the preceding section, 
Tcg/ J = 0.160±0.005. If we combine the present estimate of (3cg with the estimate 
of 7cg m the preceding section and use the scaling relations, various chiral-glass 
exponents can be estimated to be a ~ —1.7, /3cg — 1-1, 7cg — 1-5, ^cg — 1-2 
and ?7cg — 0.8. The dynamical exponent is estimated to be zcg — 4.7 by using 
the estimated value of A and the scaling relation A = /Jcg/^cg^cg- While the 
dynamical exponent zcg comes rather close to the z of the 3D Ising EA model, 
The obtained static exponents differ significantly ioxnthoiie of the 3D Ising EA 
model (3 ~ 0.6, 7 ~ 4, v ~ 2 and r\ ~ — 0.35,c3l3i2J : c3l3 suggesting that the 
universality class of the chiral-glass transition of the 3D Heisenberg spin glass differs 
from that of the standard 3D Ising spin glass. According to the chirality mechanism, 
the criticality of real spin-glass transitions is derived from that of the chiral-glass 
transition of an isotropic Heisenberg spin glass, so long as the magnitude of random 
anisotropy is not too strong. If one tentatively accepts this scenario, the present 
result opens up an interesting possibility that the universality class of many of 
real spin-glass transitions might differ from that of the standard Ising spin glass, 
contrary to common belief. 




5. Conclusion 
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In summary, spin-glass and chiral-glass orderings of isotropic 3D Heisenberg spin 
glasses are studied by extensive MC simulations. Clear evidence of the occurrence 
of a finite-temperature chiral-glass transition without the conventional spin-glass 
order is presented both by equilibrium and off-equilibrium simulations. Spin and 
chirality show very different dynamical behaviors consistent with the "spin-chirality 
separation" . While the spin autocorrelation exhibits only an interrupted aging, the 
chirality autocorrelation persists to exhibit a pronounced aging effect reminiscent of 
the one observed in the mean-field model. The universality class of the chiral-glass 
transition is different from that of the the standard Ising spin glass, while the chiral- 
glass ordered state appears to exhibit a feature of "one-step" replica-symmetry 
breaking. We expect that these numerical findings have important implications to 
the understanding of the nature of real spin-glass ordering. 
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